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Abstract 



O 

O, 

\^ i We consider an integrable scalar partial differential equation (PDE) that is second 

CN I order in time. By rewriting it as a system and applying the Wahlquist-Estabrook 

prolongation algebra method, we obtain the zero curvature representation of the equa- 
tion, which leads to a Lax representation in terms of an energy-dependent Schrodinger 
spectral problem of the type studied by Antonowicz and Fordy. The solutions of this 
PDE system, and of its associated hierarchy of commuting flows, display weak Painleve 
behaviour, i.e. they have algebraic branching. By considering the travelling wave so- 
lutions of the next flow in the hierarchy, we find an integrable perturbation of the case 
^SJ ' (ii) Henon-Heiles system which has the weak Painleve property. We perform separation 

. of variables for this generalized Henon-Heiles system, and describe the corresponding 

! solutions of the PDE. 

o 

^ - 

§ 1 Introduction 

o . 

. Recently one of us has been classifying integrable 1+1-dimensional scalar partial differential 
^ ! equations (PDEs) that are second order in time [25, including those of the type 
P. ' 

>'■ Wu = Wnx,t + F[w,Wt,W^,W^t,---] (1-1) 

X • 

^ ■ (where F is a nonlinear function of w and its x and t derivatives). This large class of scalar 
. . PDEs includes the Boussinesq equation 

Utt = + {u'^)xx, 

a well known integrable equation in shallow water wave theory. In the course of performing 
the classification, the equation 

wtt = W3^,t + 8w^w^t + 4:WxxWt - '2w^Wi^ - 4wa,j.W3a, - 24:wlw^^. (1.2) 

was found. The above equation was obtained by means of the perturbative symmetry ap- 
proach to the classification of integrable PDEs, as presented in In the latter approach, 
the existence of infinitely many commuting symmetries is taken as the defining property of 
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an integrable PDE, and by representing a PDE such as (jl.ll) symbolically (in Fourier space) 
one obtains an infinite sequence of necessary conditions for integrability. The equation ()1.2|1 
was isolated by the requirement that it should satisfy the first few of these conditions. 

Obviously, checking a finite number of necessary conditions does not prove integrability. 
However, in practice it can be proved a posteriori that only the first few conditions are 
sufficient for integrability within a given class of equations. In any case, having obtained 
a particular equation such as ()1.2|1 in the course of classifying the general class of PDEs 
(jl.lj) . its integrability must then be proved constructively, either by demonstrating that it 
is explicitly linearizable, or by showing that it can be derived as the compatibility condition 
for a linear system (Lax pair) of the form 

£^ = 0, = ^V', (1-3) 

for suitable operators C, A having coefficients dependent on w and its derivatives, and on 
a spectral parameter A. Once a Lax pair (or, alternatively, a linearization) is known, then 
the infinite hierarchy of symmetries of an integrable PDE can be obtained directly. To look 
for a Lax pair, it is most convenient to try to obtain the PDE (jl.2j) in the form of a zero 
curvature equation 

Fi-G,+ [F,G] = 0, (1.4) 
this being the compatibility condition for the matrix linear system 

= F^, ^t = G^ (1.5) 

(where the matrices F, G are functions of u and its derivatives, and of a spectral parameter). 

In the next section we apply the Wahlquist-Estabrook prolongation algebra method to 
the PDE fll.2j) in order to obtain a zero curvature representation for it. The calculations are 
much facilitated by setting u = Wx and v = Wt — w^x — and rewriting the scalar equation 
p.2j) as a two-component system (see ()2.H) below). We find that the simplest zero curvature 
representation for the system ()2.H) has F, G G sl(2), and this is equivalent to a scalar Lax 
pair ()1.3p with the Lax operator being the energy-dependent Schrodinger operator 

C = \dl+(^- + u\-\'), (1.6) 

which is an example of the general type 

N 

C = Y,^={e,dl + u,) (1.7) 

i=o 

studied by Antonowicz and Fordy [7| |H]. Some particular examples of energy-dependent 
Schrodinger hierarchies were considered earlier by Jaulent and Miodek [2H] and Martinez 
Alonso [20] (see also the more recent work of Shabat and Martinez Alonso |2I])- The Lax 
pair for p.2j) allows us to write down its hierarchy of higher symmetries. Subsequent sections 
concern the weak Painleve property for the PDE (section 3), and travelling wave solutions 
of members of the hierarchy (section 4). It turns out that the stationary flow of the next 
member of the hierarchy is equivalent to an integrable perturbation of the case (ii) Henon- 
Heiles system; this extends a result of Fordy concerning the fifth order KdV equation ^H] • 
We briefly remark on similarity reductions in our conclusions (section 5). 



2 Prolongation algebra and integrable hierarchy 

A Lax pair for the scalar PDE ()1.2|1 can be found in a straighforward way by first writing it 
in tlie form of a two- component system, that is 

Ut = U^xx + QuUx + V^, 

Vt = 4:UxV + 2uVx, 

in which form the reduction = to the KdV equation 

Ut = Uxxx + Quux (2.2) 

is made manifest. Starting from the above system, we then proceed to apply the prolon- 
gation algebra method due to Wahlquist and Estabrook jlHI (see also ^Hj for an excellent 
exposition). 

We start by assuming that the system ()2.1|) has a zero curvature representation ()1.4|) 
where the matrix F depends only on u, v, so F = F{u, v), and G depends on v and derivatives 
of u of order at most two, so G = G{u,Ux,Uxx,v) (with the dependence on a spectral 
parameter being suppressed). Plugging this ansatz into the zero curvature equation and 
substituting for Ut and Vt from ()2.1|) yields 

2 

F„ (msx + 6mMx + v^) + (Au^v + 2uVx) -G^v^-"^ G„^^ + [F, G] = (2.3) 

j=0 

(where subscripts denote partial derivatives). Reading off the coefficient of u^x in ()2.3j) . we 
have immediately that 

G = F 

SO that upon integrating we find 

G( It, Uxxy U J -f u Uxx + G{u,Ux,v) (2.4) 

for some matrix function G. Substituting back for G from ()2.4|) into ()2.3|) . we can then com- 
pare coefficients of Uxx, and we can continue in this way until we obtain explicit expressions 
for F and G as differential polynomials in u, v and their derivatives with matrix coefficients. 
These constant coefficients must further satisfy certain algebraic relations, including com- 
mutators, and for the Lax pair to exist it is necessary that there should be a non-trivial 
realization of these relations within a suitable Lie algebra. With the substitution of ()2.4|) 
and subsequent integrations and back substitutions into ()2.3|) . we find that F takes the form 

F = Au^ + Bu + Cv + D, 

where the matrices A, B, C, D are constants. However, we note that from the relations 
between A and the other coefficients of F and G, we can consistently set A = 0, so that F 
is linear in u and v. Having made this simplification, we find that G has the form 

G = Buxx + [D, B]ux + (B + 2uC)v + (6B - [B, [B, D]])y + [D, [D, B]]u + E. 

The relations that must be satisfied by the coefficients are as follows: 

[B,C]=0, 2C + [C, [D,B]] = 0, (adB)3D = 0; (2.5) 



(2.1) 



[C, (adB)2D] = 0, [B, (adD)2B] + ^[D,6B- (adB)^D] = 0; 
2[D,C] + [C, (adD)^B] = 0, [B, E] + (adD)^B = 0; 



(2.6) 
(2.7) 



[D,B] + [C,E] = 0, [D,E]=0. (2.8) 

We now seek to realize the above relations non-trivially in a suitable Lie algebra. To 
begin with we suppose that 

C = ^B (2.9) 

for some scalar ^, so that the first relation ()2.5|1 is satisfied, while the second one implies 
that 

[K,B]=2B, (2.10) 

where we have introduced 

K=[D,B]. (2.11) 

The relation ()2.1()|1 already suggests that B and K should be two of the three basis elements 
of an sl(2) subalgebra. 

If we now consider the third relation in ()2.5p then we see that 

(adB)^D = -(adB)2K = [B, 2B] = 0, 

using ()2.10|) . so this relation is automatically satisfied given the assumption (j2.9p . Further- 
more, using ()2.9p and ()2.10j) we see that both relations ()2.6p and the first relation ()2.7|) also 
hold automatically. Let us now consider the first equations in ()2.8p . which becomes 

[B,eE-D] = 0. 

The above obviously holds if we choose to set 

E = r'D, 

in which case the second relation ()2.8|) is trivially true. Finally, we have the second relation 
()2.7p still to be satified, which becomes 

-r'K+[D,[D,K]] = 0. (2.12) 

Adding a multiple of B to D does not change the commutation relation ()2.1H) . so we set 

D = -J - AB 

for some constant A, and then we find it is consistent to take J and B as the Chevalley 
generators of an sl(2) algebra, so that 

[B, J] = K, [K, B] = 2B, [K, J] = -2J, 

and the relation (|2.12|) is satisfied provided that the constants ^ and A are related according 
to 

With the above choices we see that the final expressions for F and G take the form 

'' = ("-^ + ii)«-' 



and 

G = [u^^ + 2u'^ + V + 2u\ - 4A^ + — j B + u^K - {2u + 4A) J. 
If we then choose the standard matrix representation of sl(2) then F is just 

and then setting \1/ = we see that the x part of the matrix hnear system p.Sj) is 

equivalent to the energy-dependent Schrodinger equation 

^C^^{dl + ^+u-X)i^ = (2.13) 

for ijj, with the operator C as in ()1.6p . Similarly, the t part of ()1.5|) implies that the time 
derivative of ip is given by 

dtiJ = i2u + AX)ij^-u^ij. (2.14) 

Very recently, we learnt that Shabat also derived the system ()2.H) as a reduction of the so 
called universal solitonic hierarchy; see the unnumbered equation half way down p. 622 in 

The construction of integrable hierarchies associated with energy-dependent Schrodinger 
operators of the general type (jl.7|) was described in detail by Antonowicz and Fordy [HI 13 E] • 
We now present the explicit form of their construction for the particular linear problem 
(j2.13|) . showing how this leads to the hierarchy of commuting flows for the system !^2.1\i . 
that correspond to the symmetries of the scalar PDF (jl.2j) . Starting from the Schrodinger 
equation (|2.13p . the sequence of compatible linear evolution equations for the wave function 
ip takes the form 

dt,^^, ij = {2PNd, - Pn,x) tP = Atlj (2.15) 

for integers = 0,1,2,... that label the flows, where Pn are certain polynomials in A 
whose coefficients are functions of u, v and their derivatives. For each N, the requirement of 
compatibility between the time evolution (|2.15p and the Schrodinger equation (|2.13p leads 
to the generalized Lax equation 

dt,^^,C=[A,C]+^PN,.-jC, (2.16) 
which means that the potential of the Schrodinger operator C must evolve according to 

•9*2^+1 = (dl + ^Wd, + 2W,) Pm, W = ^ + u-X. (2.17) 

For each A^, the equation ()2.17|1 encodes the evolution of u and v with respect to the time 
t2N+i, but the precise from of each P^ must be specified in order to obtain the explicit form 
of these evolution equations. In fact a compatible sequence of polynomials Pn, of degree A^ 
in A for each integer A^ > 0, can be consistently generated by considering the product 

V = tptp^ 

of two independent solutions of the Schrodinger equation (j2.13p . It is well known |25 that 
V satisfies 

VV^^- -V^ + 2WV^ + 8X = 0, (2.18) 

which is known as the Frmakov- Finney equation '■, note that we have fixed the Wron- 

skian 

i/j^ip^ - ^ptpl = aVx. 



It is also useful to observe that the differential consequence of ()2.18j) is a third order linear 
equation for "P, namely 

{dl + AWd^ + 2W^) V = (2. 19) 

We can expand the solution of ()2.18j) about A = cxo, so that 

oo 

^ = E^J-^"'' ^0 = 2' (2-20) 

j=0 

and, having fixed the sign of the leading coefficient /o, we find that all subsequent coefficients 
fj are determined uniquely from the Ermakov-Pinney equation by a recursion of the form 
fj = ( differential polynomial in M, /fc, k < j) = ( differential polynomial in m, w ). Then V 
serves as a generating function for the sequence of polynomials P/v in the sense that 



N 

N 



(A^P)+ = 5^/,A^-^- (2.21) 

j=0 



The first few of the fj are 

/o = 2, /i = n, /2 = {u,, + 3u^ + v)/4. 

In order to write down the hierarchy of flows (j2.17p explicitly in Hamiltonian form, it is 
helpful to observe that from (j2.19|) the coefficients of V satisfy the recursion 

Adjj+i = (dl + Aud, + 2u^) /, + ^ (2vd^ + V,') 

for j = 0, 1, 2, . . .. Thus, upon expanding out in powers of A, all but the two lowest orders 
in ()2.17|) cancel out, leaving the equations 



where 



and 



SuHn = = In, (2.24) 

with e.g. 

no = 2 J udx, ni = ^ J{u^ + v) dx, n2 = ^ J {-ul + 2u^ + 2uv) dx. 

These conserved quantities are just two-field generalizations of the well known Hamiltonians 
for the KdV hierarchy, to which they reduce upon setting v = 0. By the results of Antonowicz 
and Fordy concerning energy-dependent Schrodinger hierarchies [HI HE]; Bq and Bi form a 
compatible pair of Hamiltonian operators, and the existence of the sequence of functionals 
Hn satisfying (j2.24p is also guaranteed (see e.g. the Lemma on p. 468 of jS]). 

Thus we see that (|2.22j) constitutes a bi-Hamiltonian formulation of the integrable hier- 
archy whose first non-trivial member is the system ()2.ip . From fl2.24|) it is clear that the 
differential polynomials fj are gradients of conserved functionals, and the system ()2.1|) just 
corresponds to the = 1 flow with ^3 = t. We shall return to this formulation of the flows 
in section 4 when we consider travelling wave solutions. 



3 Weak Painleve property 



The connection between the singularity structure of differential equations and their integra- 
bility has been exploited since at least the time of Kowalevksi (221 • Some time later, this 
connection led Ablowitz, Ramani and Segur to make their conjecture jS] that (up to possible 
changes of variables) all ordinary differential equations arising as reductions of integrable 
PDEs should have the Painleve property, i.e. that all movable singularities of their general 
solution should be poles. While testing for the Painleve property is a very useful tool for 
isolating integrable ODEs and PDEs, it has long been known that insisting on only pole 
singularities is too strong a requirement for integrability. Indeed, there are many examples 
of integrable systems which have the weaker property that all movable singularities of the 
general solution have only a finite number of branches, and this weak Painleve property was 
proposed by Ramani et al. jHS] as a more general criterion for integrability. 

In finite-dimensional classical Hamiltonian mechanics there are many examples of Li- 
ouville integrable systems with the weak Painleve property (see e.g. Pi 121); including the 
classical problem of geodesic motion on ellipsoids that was solved by Jacobi [2Z1- In this 
section we apply the weak Painleve test of [201 to the PDE system (j2.H) . before making 
some general remarks concerning weak Painleve expansions for integrable hierarchies de- 
rived from energy- dependent Schrodinger operators, some of which were considered in [22] • 
The ODE reductions of the system ()2.1|) that are considered in the next section provide 
further examples of integrable finite-dimensional systems with the weak Painleve property. 

For the system ()2.1j) there are two types of possible expansion around movable singular- 
ities. The first type of expansion around a movable (non-characteristic) singular manifold 
(j){x, t) = has algebraic branching, and the leading order terms in this principal balance 
have the form 

« ~ 0,/(20.) + a02/^ (3.1) 

where a = a{x, t) is an arbitrary function. To find the resonances, where new arbitrary 
functions appear in the expansion, we substitute 

into (j2.H) and consider the leading order linear terms in ei and 62. This gives a 2 x 2 
homogeneous linear system for the ej, and the vanishing of the determinant of the associated 
matrix yields the polynomial equation 

_ g^3 _ ^Q^2 ^ 3^ ^ 0^ 

whose roots are the resonance values 

r = -1, 0, 2/3, 4/3. (3.2) 

The standard resonance r = — 1 corresponds to the arbitrary choice of 0, while r = 
corresponds to a as in ()3.1|) . From the fact that the other two resonances are also non- 
negative, it follows that this expansion corresponds to a principal balance, and the fact that 
these other resonances are non-integer rational numbers means that the algebraic branching 
cannot be removed simply by changes of the dependent variables (e.g. choosing as a new 
variable does not remove the cube root branching from the solution). 

The resonance conditions at the values in ()3.2|) are all satisfied, which means that the 
leading order terms (j3.1|) can be consistently extended to yield the Puiseux series 

cx> 00 

u = 0i/(20j + ^ 02/3+., ^ = J2 Vj 0-'/'+^ (3.3) 

j=0 j=0 



where the singular manifold cf) and the uq = 9vo/2 = a, U2, are arbitrarily chosen functions 
of X and t, and all other coefficients uj and vj are determined uniquely in terms of these four 
functions and their derivatives. 

The second type of expansion of the solutions u, v of ()2.H) around movable singularities 
on t) = has the leading order behaviour 

u ~ V ~ (3.4) 

where b = b{x, t) is arbitrary. In the degenerate case 6 = (corresponding to f = 0) this 
just reduces to the well known expansion for the KdV equation that was first obtained by 
the WTC method HH. This balance has resonance values 



r = -1,0,4,6, (3.5) 

where r = — 1 is standard, r = 4 and r = 6 correspond to arbitrary functions that can 
be introduced as the coefficients of 0^ and 0^ in the expansion of u (just as in the case of 
KdV). The extra resonance value r = corresponds to the arbitrary function b appearing 
at leading order in the expansion of v, as in ()3.4j) . and thus we see that this is a second 
principal balance. Also, the resonance conditions for this expansion are trivially satisfied, 
and it leads to a standard series of WTC type, namely 

oo oo 

u = 2(log 0),, + ^ = (3.6) 

i=o j=o 

where and the coefficients vq = b, U2, may be chosen arbitrarily, with all other coefficients 
being uniquely determined by them. 

It is known that for some integrable systems, there are transformations of hodograph (or 
reciprocal) type which can be used to transform a PDE with algebraic branching into one 
with the strong Painleve property P^ll5j. but it is by no means clear that this can always be 
done. In fact, in the paper |2S|; one of us made the assertion that reciprocal transformations 
of a certain type should succeed in removing all algebraic branching from the solutions of 
certain energy- dependent Schrodinger hierarchies based on Lax operators of the form 

N 

C = dl + Y,^^u,. (3.7) 
i=o 

This assertion was verified for the principal balances in such hierarchies, but a more careful 
consideration of non-principal balances suggests that the transformations considered in [231 
were actually insufficient to remove all branching from the solutions. Thus the question as to 
whether PDEs such as ()2.H) and the systems treated in [221 admit reciprocal transformations 
which transform them into PDEs with the strong Painleve property remains open. In the 
next section we consider ODE reductions of (j2.ip whose general solutions share the same 
algebraic branching as the original PDE system. 



4 Travelling waves and perturbed Henon-Heiles 



Here we consider the stationary reductions of the higher order flows of the hierarchy associ- 
ated with the system of PDEs (|2.ip . In fact, more generally we will consider travelling wave 



solutions of the higher flows ^tjjv+i such that the dependence on x and t2N+i corresponds to 
waves moving with speed c, so that 

U = u{z), V = v{z), Z = X — Ct2N+l 

(where the dependence on the other times t2j+i has been suppressed). In that case, the 
travelling wave reduction of the t2Af+i flow takes the form 



(a3 + 2(ua, + 9,n))(/^ + c/2) + i(t;9, + 9,^)/^_i = 0, 

2(vd, + d,v]ih + c/2) = 0. 



The stationary flow corresponds to setting c = 0. 



(4.1) 
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Figure 1: The profile of u{z) for a generic one-phase solution of the system h2.1\) . 

Letting f = fjsi + c/2 and g = fN-i, the second equation (|4.1|) can be integrated as 

K 

for some constant -R', and then substituting for v in the first of these equations and integrating 
again yields the differential equation 



ff"-lifr+2uf+^+L = o, 



(4.2) 



where L is another constant. 

First we consider the case = 1, which corresponds to the travelling wave reduction of 

the original PDE system ()2.1|) . In that case we have fo = g = 2, fi = u so f = u + c/2, and 
the ODE dOl) becomes 



//" - ^(/')' + 2(/ - c/2)/2 + ^ + L = 0, 
which can be further integrated to yield the quadrature 

r.M+c/2 y^J^J 



Z- Zo 



(4.3) 



(4.4) 



where Q is the quartic polynomial 



Q{f) = f- cf - Mf -Lf- K = (/ - ei)(/ - e,){f - e,){f - e,), (4.5) 

and M, zq are two more integration constants. In the affine coordinates (/, y), it can be seen 
directly that the curve 

fy' + 2Q{f) = (4.6) 

defines a two-sheeted cover of the Riemann sphere with six branch points, and so has genus 
two. This can also be seen by making the change of variables 



p2 



X = -l/f, Y = y/r 
so that the curve takes the more standard hypereUiptic form 

= 2X(eiX + l)(e2X + l){esX + l)(eiX + 1), 
and the quadrature ()4.4|1 becomes the hypereUiptic integral 

^-^o = y (4.7) 

Note that, because the inversion problem ()4.7p involves a single integral of one holomor- 
phic differential on a curve of genus two, the simple ODE (j4.3|) (which can be rewritten 
in Hamiltonian form) has deficiency one in the terminology of Abenda and Fedorov ^ |2] . 
This can be seen as the origin of the weak Kowalevski-Painleve property for the solutions. 
Nonetheless, it is possible to obtain real solutions without any singularities on the real axis, 
by choosing suitable initial data to ensure that the quantity / lies between two real branch 
points, 63 < / < 62 say. To produce numerical plots of such solutions, rather than directly 
evaluating the integral ()4.4|1 and then inverting, we have found it very convenient to do a 
numerical integration of the first order equation 



/' = ±1 



-2(/-ei)(/-e2)(/-e3)(/-e4) 



/ 



taking into account the sign changes that occur when / reaches one of the branch points 
62,63. We have plotted the profile of u{z) for the particular choice 61 = —2, 62 = —3, 
63 = —4, 64 = —5 (see figure^; this corresponds to travelling waves moving with speed 
c = Cj = —14. We refer to the travelling wave solutions given by ()4.4|1 as "one-phase" 
solutions, by analogy with the corresponding solutions of KdV and other soliton equations. 
For the system ()2.1|) it is not clear whether some such solutions have a clear interpretation 
in terms of the periodic spectral problem for the associated energy-dependent Lax operator 
()1.6|) . so the common terminology "one-gap" is perhaps not appropriate here. 

In the case when two of the roots of the polynomial Q{f) coincide, the associated hyper- 
eUiptic curve acquires a singularity and the genus drops to one; then the quadrature (j4.4p 
can be evaluated in terms of elliptic integrals. Furthermore, if two pairs of roots coincide 
(say 63 = 61, 64 = 62) then the genus of the curve drops to zero, and then the quadrature 
()4.4j) can be written explicitly as 



V2ikl - kl) ' U2 + Vkf + kl-u ) ' Wkl + kl-u + k^ 



(4. 
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Figure 2: The 'profile of u{z) for a one-soliton solution of the system \2.1\) . 



where we set 



"^2 



62 < ei = -kl < 0. 



The above choice of parameters and branches for the logarithms ensures that there is a real 
solution u = f — c/2 with 62 < / < ci < 0, so that singularities of the type f C {z — zq)^^^ 
or / ~ C (2 — zq)~^ are excluded from the real z axis. This gives a one-soliton solution for 
the system ()2.H1 . which has a kink shape (like solitons in the sine-Gordon equation). We 
have plotted an example of such a solution for the choice ci = —1, 62 = —2, with speed 
c = —6 (figure El). (Here the term "one-soliton" is used in a general sense, because as far 
as we are aware the solution does not have a direct interpretation in terms of an inverse 
scattering transform.) 

In the case N = 2 (the t^ flow of the PDE hierarchy) we find that the travelling waves 
are equivalent to the Hamiltonian system with two degrees of freedom defined by the natural 
Hamiltonian 

1 , , 9, o 1 9 £^ 256eoi , , , 

hi = ^ipl+ pI) + ql + ^qiql + 2cq, - 7^ + — j^. (4.9) 



2ql 



ql 



2''' ' ^' 2 

The coordinates and momenta are canonically conjugate, so {pjiqk\ = Sjk, and Hamil- 
ton's equations 

r\n ■ rln ■ 

{h,p,}, J = 1,2 (4.10) 



dqj _ „ I ^Pj 



dz dz 

are equivalent to the ordinary differential equation ()4.2|) for travelling wave solutions of the 
flow upon setting 



u 



9 = qi, 



4e 
P 



256e 

qt 



f = -qll^ 



4'^ 



3n' 



2c 



(4.11) 



The derivation of the equations for g^, pj in Hamiltonian form is omitted, since it is essentially 
the same as the corresponding calculation for the stationary flow of the fifth order KdV 
equation carried out by Fordy [THj; Fordy's results are reproduced by setting e = 0, c = 0. 
The equations of motion ()4.1U|) can be written in the form of a Lax equation 



where the Lax matrix L is 



with 



S(A; = 4A=^ - 2giA^ + Qg^ + Qi + A + _ ^ ^ 

The Lax equation is the compatibihty condition for the hnear system 



128egi 16e 



with 



d$ / 1 \ 

L$ = 2/x$, ^ = M$, M=(^^_^^_64e qJ. (4.13) 

The spectral curve given by det (L — 2/il) = has the explicit form 

f,' - 4A^ - 2cA3 - - h,\ - ^ + ^ = 0, (4.14) 

being the second independent integral, in involution with hi i.e. {/ii,/i2} = 0. The integral 
h2 generates a second commuting flow 

^ = {/^2,g,}, ^ = {^2,P,}, J = 1,2. 

Using the relations ()4.1ip for u and v, the joint solutions qj{z,y), Pj{z,y) of the pair of 
compatible Hamiltonian systems defined by hi and /12 yield solutions of the PDE system 
(j2.1|) upon identifying the independent variables z = x — ct^, y = it = At^. 

At this point we should comment on the integrable Hamiltonian system defined by fl4.9p . 
Homogeneous Hamiltonians of Henon-Heiles type were previously classified by using the 
usual (strong) Painleve test, requiring meromorphic solutions ^2] 5 so it is not surprising 
that extensions of such potentials with the weak Painleve property do not appear to have 
been described in detail before. However, for systems with two degrees of freedom, an 
extensive search of natural Hamiltonian systems with a second invariant has been performed 
by Hietarinta (23; ^ind he considered integrable cases of such systems with homogeneous 
cubic potentials, as well as various perturbations of such systems. The potential ()4.9p with 
a perturbation of the form qi/q2 is a special case of equation (3.2.15) in [21], and can be 
obtained by making a particular choice of the functions / and g there. 

In order to reduce the perturbed Henon-Heiles Hamiltonian system to quadratures, we 
perform separation of variables, which we now describe. The separation of variables for the 
original three integrable cases of the Henon-Heiles system (with only cubic and harmonic 
terms in the potential) was obtained in (see also ^D])- The generalized integrable cases 
(i) and (iii), which Fordy had shown to be reductions of the Sawada-Kotera and Kaup- 
Kupershmidt equations respectively (THj , were only separated quite recently jl21 • As already 
mentioned, the Hamiltonian ()4.9p that we consider here is a perturbation of the integrable 
case (ii) of the Henon-Heiles system, and so the separation of variables is straightforward. A 
new set of canonically conjugate coordinates and momenta Aj, ttj can be obtained according 
to a well known ansatz (see (221) whereby the coordinates Ai,A2 are given by the zeros of 
the upper right hand entry of the Lax matrix (j4.12j) . so we set 

4A2 + 2giA - ^ql = 4(A - Ai)(A - A2) 



to find 

2 

Ai + A2 = -|, AiA2 = -||. (4.15) 

In order to find the canonically conjugate momenta, we can write down the generating 
function 

dG dG 



G = -2pi(Ai + A2) + 4p2\/-AiA2, Qj = TTj = — , j = 1, 2, 



which means that pj are expressed in terms of the new coordinates and momenta as 
Pi \ _ V-A1A2 / -A1/V-A1A2 X2/V-X1X2 \ f T^i 



P2 J 2(Ai - A2) V 1 ~1 / V ^2 



(4.16) 



Using the above formulae, we can also rewrite tti and 7r2 in terms of the original variables 
Qj, Pj, but this will not be necessary for what follows. 

Once we are equipped with the expressions ()4.15p and ()4.16|) . we can rewrite the two 
Hamiltonian functions hi and /i2 in terms of the new variables A^ , iTj. Then one can eliminate 
between these two equations to obtain two separated equations 

A^TT^ 1 1 p 

^ - 4A^ - 2cA? - -/iiA^ - -h^X, - — + ^ = 0, J = 1, 2. (4.17) 
16 ^ ^ 2 ^ 2 ^ 64 Aj ' ^ ' 

Upon comparing the equation ()4.14|) for the spectral curve with the two separated equations 
()4.17|) . then with Hj = Ajjij/ \j the pair of points (Aj, fij) for j = 1,2 he on the spectral curve. 
Now we can write down the action S{hi, Ai, A2) such that 

dS dS 



by setting 



"''dx/ "'~dhj' 



S = r '-^ + r (4.18) 



This yields the quadratures 

f^' XdX f^^ XdX f^' dX f^^ dX 

zi= +/ , Z2= — + / — , (4.19) 

with 

zi = z + constant, Z2 = y + constant. 

The inversion problem ()4.19|) associated with the hyperelliptic curve ()4.14p of genus three 
involves only two out of three independent holomorphic differentials. This can be augmented 
with the corresponding relation in terms of a third holomorphic differential. 



dx r^^ dx 

XjJ, J AyU ' 



and then the extra variable is a function of the other two, i.e. z^ = 2:3(2:1, 2:2). So this 
perturbed Henon-Heiles system is another example of an integrable system with deficiency 
one, and instead of the standard Jacobi inversion involving the threefold symmetric product 
of the curve, one has the twofold product corresponding to pairs of points (Ai,/ii), (A2,/U.2). 
The image {zi, Z2, 23) of a pair of points under the Abel map lies on a stratum in the Jacobian 
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Figure 3: The generic shape of the profile u{x) = u{x,T) for a two-phase solution of the 
system \2.1\) . at fixed time T . 



of the curve, and this leads to an explanation for the algebraic branching in the solutions 

ID- 

The solutions of the perturbed Henon-Heiles system produce two-phase solutions of the 
PDE system (12.11) . which are the analogues of the two-gap solutions for KdV. The solutions 
corresponding to higher stationary flows (and travelling waves) of the PDE hierarchy are 
included in a general ansatz for A^-phase solutions (see jHHI, for instance). One can set 

N 

p = ^^t = -Q(A_A,), 

i=i 

where ip, ip'^ are two independent solutions of the Schrodinger equation ()2.13|) whose Wron- 
skian is 

with /i = /i(A) being given by the equation for the spectral curve, 

which is of genus + 1. Then V satisfies the Ermakov- Finney equation 

VAx - \vl + 2WV^ + = 0. (4.21) 

The expressions for u and v in terms of \j (trace formulae) are 

TV 

^ = _rf/4_2^A„ v = -^^^- (4.22) 

these follow immediately from expanding 1)4.211) around A = oo and A = 0. Similarly, 
expanding around A = Aj yields the Dubrovin equations 



]\k^j{\ - ^k) 



= ^ , , , J = l,...,iV, (4.23) 
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Figure 4: The generic shape of the profile v{x) 
system H2.1]} . at fixed time T. 



v{x,T) for a two-phase solution of the 



with prime denoting derivative with respect to x and /i being a square root of the right hand 
side of fjOnjl . 

Imposing the time dependence ()2.14j) upon ip and ip"^ means that V evolves with t ac- 
cording to 

A = {2u + AX)V., - 2u^V, 
and this imphes that the time evolution of the roots of the polynomial V is given by 

A.,. = - ^ n n r ' J = ^^---^N. (4.24) 

Just as in the one-phase case, the Dubrovin equations are extremely convenient for obtaining 
the solutions numerically. We plotted the profile u{x,t) of a two-phase solution at a fixed 
time t = T by doing a numerical integration of the equations ()4.23|) for N = 2 (see figure 
El); for this example we placed the branch points of the associated genus three curve at 
A = 1,2,3,4,5,6 (as well as the fixed branch points at A = 0, cxd). The corresponding 
two-phase profile of v{x,T) appears in figure IH 

It is also interesting to consider the two-soliton solutions of ()2.H) which arise from the 
coalescence of three pairs of branch points on the genus three curve. In order to produce 
plots of such a solution, we performed numerical integration of the Dubrovin equations 
(j4.23|) for = 2, in the particular case when the right hand side of (j4.2(J|) has double roots 
at A = 1,2,3, using particular initial conditions Ai|2,.=o = Ai(0,To), X2\x=o = -^2(0, To), with 

l<Ai(0,To)<2, 2<A2(0,To) <3. 

This produced the two-soliton profile of u{x,Tq) = — 2Ai(x, Tq) — 2A2(x, Tq) — d/i with 
d = —48, corresponding to the configuration of two solitons at some fixed time t = Tq 
(plotted in figure E))- In order to observe the subsequent evolution of the pair of solitons at 
times t > To, we then integrated the two equations ()4.24j) for j = 1,2, starting from the 
initial data lai\t=To = Ai(0,To), la2\t=o = Ai(0,To), in order to get values for Aj(0,t) for 
t > Tq. For a fixed sequence of times t = Tn, the values Aj(0, T„), j = 1, 2 were then used as 
initial data for the Dubrovin equations (j4.23|) . in order to produce the sequence of profiles 
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Figure 5: The profile of u{x) = u{x,Tq) for a two-soliton solution of the system at 
time Tq (before collision). 
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Figure 6: The profile of u{x) = u{x,Ti) for a two-soliton solution of the system h2.1\] at 
time Ti (just before collision). 







4.4 












/ ^ 












1 4. 


- 






-10 


-5 




C 


) 5 

X 


10 



Figure 7: The profile of u{x) = T2) for a two-soliton solution of the system \2.1\) at 
time T2 (during interaction). 

of u{x, Tn). The interaction between the two sohtons can be seen in figures IHl- El The shape 
of the field v{x,t) = e/(A^A2) for the two-sohton solution (with e = —144 in this example) is 





qualitatively similar to that of u{x,t), but for completeness we have also plotted v{x,T2) in 
figure EH which is at the same stage of the soliton interaction as the corresponding profile 
of u in figure 
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Figure 8: The profile of u{x) = T3) for a two-soliton solution of the system \2.1\) at 
time T3 (just after collision). 
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Figure 9: The profile of u{x) = u{x,T^ for a two-soliton solution of the system h2.1\l at a 
later time T4 (after interaction). 



5 Conclusions 

The Wahlquist-Estabrook method has proved to be effective for finding the Lax pair for the 
PDE p.2|) . leading to a 2 x 2 zero curvature representation. The latter is equivalent to a 
scalar Lax pair involving the energy-dependent Schrodinger equation ()2.13p . The solutions 
of integrable hierarchies associated with energy-dependent Schrodinger operators have some 
interesting features compared with those of other integrable systems jH E]- As far as we 
are aware, the theory of inverse scattering has not been developed for such operators (apart 
from the case where C is quadratic in A, starting with the work ^^). 

The solutions of such hierarchies generally display the weak Painleve property and in 
particular the algebro-geometric solutions often correspond to finite-dimensional integrable 
Hamiltonian systems with deficiency p. For the iV-phase solutions of (j2.H) described above, 
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Figure 10: The profile of v{x) = v{x,T2) for a two-soliton solution of the system \2.1\) at 
time T2 (during interaction). 

the deficiency is one: the genus of the curve ()4.20j) is one higher than the number of de- 
grees of freedom. It is interesting to note that, while many such Hamiltonian systems were 
obtained in j2] by applying Dirac reduction to algebraically completely integrable systems 
that linearize on a (jf-dimensional Jacobian associated with a genus g curve (with the number 
of constraints being the deficiency), the systems considered here can be viewed instead as 
deformations where the genus increases from g to g + 1: the parameter e is the deformation 
parameter, and the spectral curves for the KdV hierarchy are recovered in the limit e — > 0. 
In particular, the perturbed Henon-Heiles system with Hamiltonian ()4.9p corresponds to a 
deformation from genus two to genus three. One can obtain multi-parameter deformations 
(with correspondingly greater deficiency) by considering the Lax operator 



We have plotted some one- and two-phase solutions of the PDE system ()2.H) . includ- 
ing their solitonic limiting cases. The soliton dynamics has some quite different features 
compared with KdV solitons, and is worthy of further study. In the absence of an inverse 
scattering transform for the PDE system ()2.1|) . one can wonder to what extent the class 
of soliton-type or A^-phase solutions represent the general solution. For the case of some 
integrable PDEs, such as the KdV and nonlinear Schrodinger equations [201; it is possible 
to make local approximations to arbitrary periodic solutions using finite-gap ones. 

Since the first draft of this paper was written, it has been brought to our attention that 
the Lax pair for the system ()2.H1 was written down by Zakharov some time ago jlHj. Ito 
found the bi-Hamiltonian structure and conservation laws for this system rewritten in the 
variables u and V = y/v [22] • In this form, Matsuno studied the dispersionless limit and also 
gave an implicit expression for the travelling waves of the system, providing a parametric 
formula for the quadrature ()4.4j) . or equivalently the hyperelliptic integral ()4.7p . in terms 
of Jacobi elliptic functions (see Appendix B in [32] )• More recently, the system ()2.H) has 
reappeared in a symmetry classification of mixed scalar and vector evolution equations (see 
HU, section 3.2.3). 

Scaling similarity reductions of integrable PDEs commonly yield non-autonomous ODEs 
of Painleve type, and it has been observed that these often inherit the same Hamiltonian 
structures as their autonomous counterparts (see e.g. [221 and references). However, for 
scaling similarity reductions of the PDE hierarchy associated with (j2.H) this appears not to 
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be the case. For example, following the method of Clarkson and Kruskal [M] we see that by 
setting 

u{x,t) =e^U{z), v{x,t) = 9^V{z), 

with 

(so that 6 = {—3t + 3to)"^^^), (j^.ip is reduced to a coupled system of ODEs, namely 

U'" + 6UU' + V - zU' - 2U = 0, 

{2U - z)V' + A{U' - 1)V = 0, ^ ' 

where the prime denotes d/dz. The integrating factors which allow the ODEs for travelling 
waves to be written as the single second order equation ()4.3j) no longer work for this non- 
autonomous system, except in the special case = when ()5.H) reduces to the equation 
PXXXIV (see j2S], chapter XIV). Aside from this special case, the general solution of the 
fourth order system ()5.1|) will have the weak Painleve property. Nevertheless, by applying 
the technique of Flaschka and Newell ^7j, we find that putting 

0(^;C)=^V(a:,t;A), C = 0-'\ 

leads to the linear system 

0" + (f/-C + V^/(4C))</' = 0, 

(d^cj) = _(2f/-^ + 4C)0' + (f/'-2)0. ^ ' 

The fourth order ODE system arises as the compatibility condition for this linear system, 
so while (15. 1|) does not have the Painleve property it should still be possible to interpret it 
via isomonodromic deformation of the second equation ()5.2j) . 

Another integrable scalar equation of the type p.ip considered in [33] is equivalent to 
the system 

ut = U5x + '25ujcUr,x + 10uu3r, + 20u\x + V:^, (5.3) 

The above system is a rewriting of equation number (101) in [33], which is obtained by taking 
u = Wx and eliminating v. This system reduces to the Kaup-Kupershmidt equation when 
w = 0. By analogy with our results for the system (jl.2p . related to an energy- dependent 
operator of second order, it is natural to guess that the system ()5.3|) should possesses a scalar 
Lax pair with an energy- dependent Lax operator of third order, which generalizes the Lax 
pair of the Kaup-Kupershmidt equation. Indeed this turns out to be the case, and we obtain 
the Lax pair 

\cil) = (dl + 2udx + - - a) = 0, (5.4) 



V " ^ ^ 9A 
dti) = -9Xtpxx + {uxx + 4u^)V'x - (^i3x + ^uux + I2u)il) (5.5) 

for ()5.3p . This Lax pair is one member of a family classified by Antonowicz, Fordy and Liu 
p. It is also equivalent to the 3x3 zero curvature equation recently given in [2Z1- By a 
slight extension of Fordy's results in ^H], we find that the travelling wave reduction of ()5.3|) 
corresponds to another perturbed Henon-Heiles system, with Hamiltonian 

h = \{pl+ pI) + ^q! + \qiql + 2cgi - ^ + ^. (5.6) 



This reduces to the (generahzed) integrable case (iii) Henon-Heiles system when e = 0, for 
e 7^ it has the weak Painleve property. The possibihty of including the extra term l/gf 
in the potential, while still preserving integrability, was noted some time ago by Hietarinta 
|21j . However, at present we do not know how to carry out separation of variables for this 
perturbed Hamiltonian system; the results of |121 cannot be directly extended to this case. 
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